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Abstract 



We show that in the point process limit of the bulk eigenvalues of /3-ensembles of 
random matrices, the probability of having no eigenvalue in a fixed interval of size A 
is given by 



(,. , + ,.( I DA' . ,.xp (-^A 2 + 



2-i|A 
8 4 



as A — > oo, where 



7/3 



2 13 



This is a slightly corrected version of a prediction by lDysonl (|l962T ). Our proof uses the 
new Brownian carousel representation of the limit process, as well as the Cameron- 
Martin-Girsanov transformation in stochastic calculus. 



1 Introduction 



In the 1950s Wigner endeavored to set up a probabilistic model for the repulsion between 
energy levels in large atomic nuclei. His first model s were r a ndom mer omorphi c func- 
tions rel ated to random Schrodinger operators, see IWignerl (|195ll) and IWignerl (|1952l ). 
Later, in IWignerl (|1957|) he turned to models of random matrices that are by now stan- 
dard, such as the Gaussian orthogonal ensemble (GOE). In this model one fills an n x n 
matrix M with independent standard normal random variables, then symmetrizes it to 
get 



.4 



M + M T 



The Wigner semicircle law is the limit of the empirical distribution of the eigenvalues of 
the matrix A. However, Wigner 's main interest was the local behavior of the eigenvalues, 
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namely the repulsion between them. One mathematical manifestation of this repulsion is 
large gap probabilities. If we scale the eigenvalue point process of A so that the average 
spacing is 2tt, then Wigner predicted that the probability that there is no eigenvalue in a 
fixed interval of length A is given by 

p A = exp(-(c + o(l))A 2 ), 

where this is a A — > oo behavior, and we assume that the n — > oo limit has already been 
taken. This rate of decay is in sharp contrast with the exponential tail for gaps between 
Poisson points; it is one manifestation of the more organized nature of the random eigen- 
values. Wig ner" s estimate of the constant c, l/(16n), later turned out to be inaccurate. 
Dysonl |l962) improved this estimate to 



P\ = («/3 + o(l))A 7 ' 3 exp 



64 







A 



(1) 



where (3 is a new parameter introduced by noting that the joint eigenvalue density of the 
GOE is the (3 = 1 case of 



~/3£IUA 2 fe /4 



A,r 



(2) 



j<k 



The family of distributions defined by the density (O is called the /5-ensemble. Dyson's 
comput ation of the exponent ^g, namely ~(| + | + 6), was shown to be slightly incorrect. 



Indeed, Ides Cloizeaux and Mehtal (|1973l) gave more substantiated predictions that 7/3 is 
equal to— 1/8,— 1/4 and —1/8 for values (3 = 1, 2 and 4, respectively. Math ematically pre - 



Deiftetal 



rise proofs for th e (3 = 1, 2 and 4 cases were later given by several authors: IWidoml (fl996l ). 



19970 . Moreover, the va lue of Kg and hig h er order asymp t otics were a 



Deift et al. 



so es- 



(|2007b . 



tablished for these specific cases by lKrasovsky ( 2004 ), Ehrhardt ( 2006 ), 
The problem of determining the asymptotic probability of a large gap naturally arises 
in other ran dom matrix models a s well. For a treatment of the case of the /3-Laguerre 
ensemble see lChen and Manning! (|19960 . 



Our main theorem gives a mathematically rigorous version of Dyson's prediction for 
general (3 with a corrected exponent 7 p. 

Theorem 1. The formula flU holds with a positive Kp and 

The proof is based on the Brownian carousel, a geometric representation of the n — >• 00 
limit of the eigenvalue process. We first introduce the hyperbolic carousel. Let 
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• b be a path in the hyperbolic plane 

• z be a point on the boundary of the hyperbolic plane, and 

• / : R + — > M + be an integrable function. 

To these three objects, the hyperbolic carousel associates a multi-set of points on the real 
line defined via its counting function N(X) taking values inZU{— 00,00}. As time in- 
creases from to 00, the boundary point z is rotated about the center b(t) at angular speed 
Xf(t). N(\) is defined as the integer-valued total winding number of the point about the 
moving center of rotation. 

The Brownian carousel is defined as the hyperbolic carousel driven by hyperbolic 
Brownian motion b (see Figure 1). It is connected to random matrices via the following 
theorem: 



Theorem 2 (jValko and Virag ( 2008b ). Let A n denote the point process given by (0), and let jx n be 

a sequence so that n^(2y/n— \fi n \) — > 00. Then we have the following convergence in distribution: 



y/4n- /4(A n - fj, n ) =^ Sine^, 
where Sine^ is the discrete point process given by the Brownian carousel with parameters 



(3) 



/(*) 



(4) 



and arbitrary z. 





Figure 1: The Brownian carousel and the winding angle a x 
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Remark 3. The semicircle law shows that most points in A n are in the interval [—2y/n, 2y / n]. 
The discrete point process A n has two kind of point process limits, one near the edges of 
this interval and another in the bulk. The condition on the parameter \i n means that we get 
a bulk-type scaling limit of A n . The scaling factor in © is the natural choice in view of the 
Wigner semicircle law in order to get a point process with average density l/(2n). The 
limiting point process for the edge-scaling case have been obtained by Ramirez, Rider, 



and Virag (|2007l) 



The Brownian carousel description gives a simple way to analyze the limiting point 
process. The hyperbolic angle of the rotating boundary point as measured from b(t) fol- 
lows the following coupled one-parameter family of stochastic differential equations 

da x = Xfdt + Re((e~ iax - l)dZ), a A (0) = 0, (5) 

driven by a two-dimensional standard Brownian motion and / given in ©. For a single 
A, this reduces to the one-dimensional stochastic differential equation 

da x = Xfdt + 2 sin{a x /2)dW, a A (0) = 0, (6) 

which converges as t — ► oo to an integer multiple a\(oo) of 2n. In particular, the number 
of points of the point process Sine/? in [0, A] has the same distribut ion as 0^(00) /(27r) and 



Px is equal to the probability that a converges to as t — > 00. See lValko and Viragl (|2008l) 
for further details. 

In the analysis of equation ([6]) it helps to remove the space dependence from the diffu- 
sion coefficient by a change of variables X(t) = log(tan(a(i)/4)). The diffusion X satisfies 
the stochastic differential equation: 

dX = ^fcoshXdt + 1-tanhXdt + dB, X(0) = -00. (7) 

_ _ 



In 



Valko and Virag ( 2008^ equations © and were used to identify the leading term in 



the asymptotic expansion of p x in ([!]). The proof of Theorem [T] requires a more careful 
analysis of equation ©. 

In Lemma H] we will show that for any initial condition X(0) = x £ [— 00, 00) there is a 
unique solution of the equation given in © and p x = p\(— 00) where 

Px(x) := P (X(t) is finite for all i > and does not converge to 00 as t — >• 00) . (8) 

A time shift of equation © only changes the parameter A and the initial condition. This, 
together with the Markov property of the diffusion X, shows that with T — | log A we 
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have 



Px = E 



1 {X(t) is finite f or all < t < T } ■ pi(X(T)) 



(9) 



Our main tool is the Cameron-Martin-Girsanov formula, which allows one to compare 
the measure on paths given by two diffusions. If we knew the conditional distribution of 
the diffusion X under the event that it does not blow up, then we could use the Cameron- 
Martin-Girsanov formula to compute p\ explicitly. While we cannot do this, the next 
best option is to find a new diffusion Y which approximates this conditional distribution. 
The density (i.e. the Radon-Nikodym derivative) of the path measures given by Y with 
respect to the measure given by X will be close to the right hand side of dU). Our strategy 
for finding Y is described in Section HI 

In Section|5l we will present a coupling of the transformed processes that enables us to 
show that the asymptotics is precise up to and including the constant term Kp. The term 
Kp is then identified as the expectation of a functional of a certain limiting diffusion. 

Open problem 1. Give an explicit expression for Kp for general values of 3. 

The known values of Kp are 

Kl = 2 13/24 e |C'(-l) ) K2 = 2 r/12 e 3 C '(-l) ) ^ = 2 -13/12 e |C'(-l) ) 

where £'(— 1) is the coefficient of the linear term in the Laurent series of the Riemann-£ 
function at -1. 

A natural generalization of Theorem[T]would be to consider the asymptotic probability 
that there are exactly n eigenvalues in a large interval [0, A]. This probability is usually 
denoted by Ep(n : X) in the li t eratu re. For 3 = 1, 2 and 4 the following large A asymptotics 



was obtained by lBasor et al.l (|1992l) : 



Tio ti l 3 ti3 
log Ep(n; A) = log^(0; A) + -fX + - f 1 - | - -f ) log A + Cp + o(l) 



with explicit constants cp. (See also iTracy and Widoml (|1993h .) We believe that our meth- 



ods can be used to extend the previous asymptotics for general values of 3. 

The rest of the paper is organized as follows. In the next section, we justify the p\ = 
px(-oo) and give some preliminary estimates on the probability p\(r) appearing in ©. 
Section [3] presents the version of the Cameron-Martin-Girsanov formula that we need. 
Section H] describes the strategy for finding Y and Section [5] builds on these sections to 
complete the proof of the main theorem. 
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2 Preliminary results 



First we formally verify the connection between the gap probability p\ and the diffusion 
given in (0. 

Lemma 4. The diffusion © has a unique solution for any initial condition X(0) =16 [—oo, oo) 
and p\ = px(-oo). 

Proof. The change of variables function logtan(-/4) is one-to-one on (0, 2%) — > M. There- 
fore, even with — oo initial condition, the diffusion X is well defined and has a unique 
solution until a reaches 2ix, when it blows up. We define X{t) = oo after this blowup. 
Note that for A > the solu tion of equation© is always monotone increasing at 



multiples of 2n. See Section 2.2 in lValko and Viragl (|2008h for more details. So if a(t) — ► 



as t — > oo then < a(t) < 2n for all t > 0. This means that X{t) is finite for all t > and 
X{t) cannot converge to oo which proves p\ = p\(—oo). □ 

Next, we prove a preliminary estimate on the blowup probability of the diffusion (0. 

Lemma 5. Recall that pi(x) is the probability that the diffusion with A = 1 and initial con- 
dition X(0) = x does not blow up infinite time and does not converge to +oo as t — > oo. We 
have 

< px{x) < c^exp 

Proof. For the upper bound we first assume that x > 4. Consider the diffusion 

dR= —e R ~ mt dt + dB, R(0)=x. (10) 
16 

This has the same noise term as X. The drift term of R is less than the drift term of X 
(with A = 1) as functions depending on time and space when the space variables are 
nonnegative. Thus while R stays positive, we have R < X. This means that for every 
t > we have 

Pi(x) < P (X does not blow up before time t) 



< P min R(s) < or R does not blow up before time t . (11) 

\se[o,t] * 



The difference Z = R — B satisfies the ODE 





e- z dZ = ^ e B - p/it dt, Z(0) = x. (12) 
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Integration gives 

e -x_ e -Z{t) = R_ f B{.)-fi/A. d 

16 Jo 

This shows that Z is increasing in t, in particular Z{t) > x. So if minr 0) t] R < then 

mini? < —x. 

[o,t] 

Furthermore, if 

e~ x <^- FeBM-W'ds (13) 
16 Jo 

then R blows up before time t. This certainly happens if the minimum of B on the interval 
[0, t] is not sufficiently small. More precisely if 

£_(! _ e W4) > e -« (14) 

and miri[ ,t] B > —b then (fl3~|) follows. So if 6 < x, and (p!4|) holds, then the right hand side 
of ((IT) can be bounded above by 

P fminfi < -6^ = P(|B(t)| > b) < ^ . 
\[o,t] J b 



. 16 2-x z. 4e 



We set 

/3 " ' " v/30 
As x > 4, both 6 < a; and (fi~4]) are satisfied and we get the upper bound 

Pi{x) < — e ™<cpe so e 

with C/ 9 = The upper bound for all values of x now follows by changing the 

constant cp appropriately. 

For the lower bound note that since the Sine/? process is discrete and translation in- 
variant in distribution, there exists v G (0, 1) so that p v = p v (—oo) > 0. By the Markov 
property, we have 

/CO 
K ,i(-oo, dx)p ve p/4(x) 
-oo 

where K Stt (y, dx) is the transition kernel of the Markov process X with parameter A = v. 
This implies that for some i Gl we have 

Consider the process X started at x with parameter A = 1. The Markov property applied 
at time t = 1 — | log v and the monotonicity of p\(x) in x implies 

Pi(x) > P (X(t ) < x ) p ue p/i(x ) > 0, 

since P(X(t) < x) is positive for all x e R and t > 0. □ 
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3 The Cameron-Martin-Girsanov formula 



The diffusions we work with can blow up to +00 in finite time, in which case they are 
required to stay there forever after. For this reason, we will not be able to use the text- 
book versions of the Cameron-Martin-Girsanov formula, and have to prove a version that 
applies to such diffusions. 

Proposition 6. Consider the following stochastic differential equations 

dX = g(t,X)dt + dB, limX(f) = -00 (15) 
dY = h(t,Y)dt + dB, limY(t) = -00 (16) 



on the interval (0,T], and assume that (|15|) has a unique solution X in law taking values in 
(—00, 00]. 
Let 

G S = G S (X)= [ h(t,X)-g(t,X)dX-± [ h(t,X) 2 -g(t,X) 2 dt (17) 
Jo ^ Jo 

and assume that 

(A) g 2 — h 2 and g — h are bounded when x is bounded above. (Then G s is almost surely well- 
defined when X s is finite.) 

(B) G s is bounded above by a deterministic constant. 

(C) G s — > —00 when s|r if X hits +00 at time r. In this case we define G s := —00 for s > r. 

Consider the process Y whose density with respect to the distribution of the process X is given 
by e° T . Then Y satisfies the second SDE (jl6|) and never blows up to +00 almost surely. Moreover, 
for any nonnegative function (p of the path of X that vanishes when X blows up we have 

E^(X) = E[^(Y>- GT(y) ]. (18) 

Proof. We follow the standard proof of the Girsanov theorem. 

First we show that G s is well defined for finite X s . From condition |(A)| it follows that 
the second integral is well-defined. The first integral can be written as 

PS 

(h-g)dB+ / (h-g)gdt 
Jo 

which is well defined since (h — g) and 2{h — g)g = [h 2 —g 2 ) — (h — g) 2 when their argument 
x is bounded above. 
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Next, we show that M s = e Gs is a bounded martingale. This is clear after the hitting 
time r of X of +00, if such time exists. Before this time, G s is a semimartingale, and so is 
M s . Ito's formula gives 

dM =(h- g)MdB 

so that the drift term of M vanishes. So M is a local martingale which is bounded, so it 
has to be a martingale. 

The rest of the proof is standard, and follows the following outline. Set 



PS PS 

B S = X S - I hdt = B s - (h- g)dt. 
Jo Jo 



It suffices to show that B is a Brownian moti on with respect to the new m easure with den- 
sity Mt- This follows from Levy's criterion (IKaratzas and Shrevel (|199ll) Theorem 3.3.16) 
if B and B 2 — s are local martingales. Since M is a martingale, it suffices to show that 
BM and (B 2 — s)M are local martingales with respect to the old measure, which is just a 
simple application of Ito's formula. 

The identity (|T8|) is just a version of the change of density formula. □ 



4 Construction of the diffusion Y 

In this section we will create a diffusion which approximates the conditional distribution 
of the diffusion X under the event that it does not blow up. We will construct a drift 
function h(t, x) for which the diffusion Y 

dY = h{t, Y)dt + dBt, Y(0) = -00 (19) 

is well-defined, a.s. finite for t > and the (formal) Radon-Nikodym derivative e° T with 
Gt defined in (|T7)) is almost equal to the right hand side of equation (d) with the appro- 



priate 7/3. 

Lemma 7. For the diffusion 0, A > 1 and T = flog A there exists a function h(t,x) so that 



conditions \(A)§(C)\ ofProposition\6\hold, and G T has the following form: 



-G T (X) = -^A 2 + (^-\) \ + l(p+^-6)\og\ (20) 

T 

<p(T-t,X(t))dt. 



64 V § 4 / 8 V P 
+^ + (2-^jX(T)+ + u(X(T)) + j o 



Here the function u is bounded and continuous, 4> is continuous and \4>(t,x)\ < 4>{t) with 
Jo°° ^(t)dt < 00. The functions u and </>, <p may depend on the parameter (3, but not on A. 
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The function h will have the following form 

h(t, x) = — — /sinh(x) + ho(t, x), (21) 
where \h (t, x)\ < cifO < t < T. The constant c depends only on f3. 

Proof. Construction of the function h. Given an explicit formula for h it would not be 
hard to check that Gt has the desired form. However, we would like to present a way 
one can find the appropriate drift function. This will provide a better understanding of 
the form of the resulting h. 
We will use the definition 

-G S (X)= [ g(t,X)-h(t,X)dX + - [ h 2 (t,X) - g 2 (t,X)dt. 
Jo 2 J 

where 

9 = 9i + 92, 9i(t,x) = -/(*)coshar, g 2 (t, x) = -t&nhx. 

Our goal is to find the appropriate drift term h in a way that the diffusion Y will approxi- 
mate the conditional distribution of X given that it does not blow up in the interval [0, T). 
We will do this term by term, starting with the highest order; towards this end we write 

h = hi + h 2 + h 3 + h 4 . We set 

h 1 {t,x) = ~f(t)smh(x) (22) 
as this yields the nice cancellation 

hi - 9\ = * f(tf sinh 2 (X) - £ f{tf cosh 2 (X) = f{tf 

in the main terms of h 2 — g 2 . In addition, if the remaining term h 2 + h s + h 4 is bounded, 
then it will be easy to show that conditions (A) (C) of Proposition [6] are satisfied. This will 
be done at the end of the proof. 

The contribution of the drift terms h% and g\ to the stochastic integral part of — G s is 
given by 

^ J° /(t)(coshp0 - smh(X))dX = ^j^ e x f(t)dX. (23) 

Our main tool for evaluating integrals with respect to dX is the following version of Ito's 
formula. Let a, b be continuously differentiable functions and let a denote the antideriva- 
tive of a. Then 

a(t)b(X) dX = d{a{t)~b{X)) - a'{t)b(X) dt - -a{t)b'{X) dt. (24) 
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Since f'{t) = —f3/Af(t), and X(0) = — oo, this formula gives 



A 



A 



A f(3 1 



- / f(t)e*dX=-f(s)e^ + -r--- 



e x fdt. 



(25) 



Next we would like to choose h 2 in (|22|) so that the integral term in the right hand side of 
((251) simplifies. More precisely, since we expect the diffusion X to be near most of the 
time, we would like to replace the term e x by 1. The plan is to use the cross term J h\h 2 dt 
in the | J h 2 dt term of G to do this. Namely, we would like to have 



h x h 2 

The solution for (|26l) is given by 



A ff3 



2 V4 



;i - e x )f- 



(26) 



//,(/.. r) • ( J - I ) ( I +tanh(.r/2)). 



(27) 



We will choose the next term, h 3 , so that the cross term J hih 3 dt in | J /i 2 eft cancels the 
cross term — J g x g 2 in — | f g 2 dt. This leads to the equation 



hh = gm = ^ f{t) cosh(x) • i tanh(x) 



which gives 



h 3 (t,x) 



(28) 



Collecting all our previous computations we get 



dt 



+ - / 2/ii/i 4 + (/^2 + + M 2 - g\ dt 
2 Jo 

— / h 4 dX + g 2 — h 2 — h 3 dX. 



(29) 



The integrand u = g 2 — h 2 — h 3 in the last integral of (|29|) has antiderivative 



/9 



■u(x) =(l x+ll log cosh(x/2) H — log cosh a;. 



(30) 



By Ito's formula, J* u(X)dX — u(X)\ s Q is given by 
u'{X)dt = 



1 

~2 
(3-6 
16 



- — - sech(X/2) 2 + - sech(X) 2 
8 2 



eft 



s + 



(2-0) ^,^2 1 



tanh(X/2f + -tanh(X)' 



(31) 
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Note that 

lim u{x) = — — -log 2 = c\. 

X— > — oo 2 

Substituting this computation for the last integral and expanding (h 2 + h 3 + h 4 ) 2 we can 
rewrite d29l) as follows. 



4/ ./n \ 2 V 4 / 16 

+ ^/(s)e*« (32) 



+ - / (2/ii/i 4 + 2(h 2 + h 3 )h 4 + hi) dt- h 4 dX + / 
2 Jo Jo Jo 



rj{X{t))dt. 



The coefficient of s in the first line of ([32]) comes from the first term on the right in ((31) 
and the constant term of (h 2 + h 3 ) 2 /2. The function r] collects the terms from the integrand 
in (|3l]), the terms (h 2 + h 3 ) 2 /2 with the constant term (/3/4 - l) 2 /2 removed, and -g 2 /2. 
More explicitly, we have 

r}(x) = ^ - 6 f n +l3 - (2 tanh(ar/2) + tanh(x/2) 2 ) + i(tanhx) 2 . 
32 8 

The function rj(x) contributes to an error term that needs to be controlled, but whose 
precise value does not influence our final result. Now we are ready to set the value for 
h 4 : we will choose it in a way that the cross term J h\h 4 dt in (|32|) will cancel the integral 
J 7] dt. This gives h 4 = —r)/hi, that is 

h 4 = 2 ^ (33) 
Xf(t) sinri(x) 

The function h 4 is a product of a function of t and a function of x. Ito's formula (f2~4~l> , with 
the notation h 4 (t, x) = f* h 4 (t, y)dy yields the evaluation of the stochastic integral in (|32|) : 

-f h 4 dX = -h 4 (s,X(s)) + j [ h 4 dt+~ f d x h 4 dt. 
Jo 4 J 2 J Q 

Plugging this into (|32|) and simplifying the deterministic terms in the first line of (|32~1) we 
arrive at 

-G, = -4- 



fdt + X^- 1 -) J* fdt + ±(f3 2 + 12P + 8)s 
+Ae-^ 4s f e x{s) + u{X(s)) - a - h 4 (s, X(s)) (34) 

o 

+ J f 2(/i 2 + h 3 )h 4 + h\ + jh 4 + ^d x h 4 ^j dt. 
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Note that h 2 and h 3 do not depend on t and are bounded by an absolute constant. The 
functions /i 4 , J14, d x h± are all bounded by a constant times l/(A/(t)) = jpf(T — t), which 
itself is bounded by a constant not depending on A as long as < t < T. Thus we can 
rewrite the integrand in (|34"1> as 

f s 

4>{T -t,X(t))dt (35) 







with a continuous function <\> which does not depend on A and satisfies \<j)(t, x)\ < 4>(t) 
with J °° (j)(t)dt < 00. Using (|3Q]> and the fact that log cosh x — \x\ is bounded, the terms in 
the second line of (|34|) can be written as 

(T)) (36) 



(2 - X(s)+ + \e~^e x ^ + wb(X( S )) - h 4 (T, X 



with a bounded and continuous ujq. This concludes the construction of the function h. In 
order to get the expression (120) for — Gt we first plug in s = T into ((34]) . Then the first line 
gives 



X 2 p 



64 

and by (|36|) the second line transforms to 

(2 - X(s) + + j?e*M + ^ (X(T)) - ^(T, X(T)). 

Note that the expression /i 4 (T, x) does not depend on T and is bounded. This proves that 
— Gt is in the desired from (|20|) . 

Now we are ready to check that the proposed choice of h satisfies all the needed con- 
ditions (A)-(C). 



Condition (A) As x — > — 00 we have 



1 _ 1 1 _ 1 

g(t, x) = -Xfe x - - + g(t, x), h(t, x) = -Xfe x - - + h(t, x) 

where \g\ < ce x and \h\ < c'e x with constants that only depend on (3 if < t < T. From 
this it follows that g — h and g 2 — h 2 are both bounded if x is bounded from above. 



Condition |(B)j We need that (|34|) is bounded from below if < s < T. The integrals 



in the first line are bounded by a constant depending on A and j3 only. The same is true 



for the integral in the last line, see (|35|) and the discussion around it. Thus we only need 



to deal with the evaluation terms of the second line. By (|36|) we just need to show that 

2 - X (s) + + Ae-^e x W (37) 
2 J 8 
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is bounded from below. Since s < T = j- log A, we get that ((371) is bounded from below by 

P\ vt _\ + , P _X(s) 



which in turn is bounded from below by a constant depending only on (5. 



Condition |(C)1 This follows the same way: one only needs to check the behavior of 
((37} as s converges to the hitting time of oo. This expression converges to oo as X(s) — ► oo 
which means that G s — > — oo. □ 

5 The proof of the main theorem 

We are ready to prove Theorem [1] 

Proof of Theorem [I] Lemma |4] gives p\ = p\ ( — oo ) , where 

p\(x) = P (X(t) is finite for all t > and does not go to oo as t — > oo) 

with X(0) = a;, as defined in ©. Note that a time shift of equation (0 only changes A and 
the initial condition. With 

T = T x = Uog\ (38) 

the diffusion r i— > X(r + T) satisfies © with A = 1 and with initial condition — oo at 
r = —T. This suggests that the dependence on A for the probability on the right hand 
side of © comes mainly from the interval [0, T\. Because of this we take conditional 
expectations in © with respect to the cr-algebra generated by (X(t),t e [0, Tj). Using the 
Markov property of X we obtain 

Px = E (1 {X(t) is finite for all < t < T } ■ Pl (X(T))) (39) 

The first term in the expectation is a function of the path X(t) on the time interval [0, T). 
Consider a diffusion Y given by the SDE (|T9|) with a drift function h(t, x) given by Lemma 
[3 With the notation of Lemma [7| we set 

rl>0O = ( 2 " f ) Y ( T ) + + ^ Y(T) + + <KT - t, Y(t)) dt. (40) 



We apply the Girsanov transformation of Proposition [6] together with equation (|20|) of 
Lemma [7] to get 

p x = A^e"^ 2 +(f-i) A E [ Pl {Y(T x )) exp{^(F)}] , 
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where 7^ = | + | — 3V In order to prove the theorem it suffices to show that the limit 

lim E[ Pl (F(T A ))exp{^(F)}] (41) 

A^oo 

exists, and is finite and positive. This limit then equals the constant Kp of the asymptotics. 
Recall that in (|40l> the function to is continuous and bounded and <f>(t, y) can be dominated 
by a function <fi(y) which has a finite integral in [0, 00). 

We will run the process Y\(t) with a shifted time, t = t — T = t — -| log A; that is, let 

Y T (r) :=Y x (r + T). 

The advantage of this shifted time is that the diffusions Y T (r) for different A satisfy the 
same SDE except they evolve on nested time intervals: 

dY T (r) = h(r, Y) dr + dB, r > — T, Y T (-T) = -00, (42) 

where the drift term is given by 

h(r, y) = h(T + T,y) = -^e"^ 4 sinh(y) + h (T + r,y). (43) 

In this new time-frame we need to show that the limit 

lim Ep 1 (F T (0))exp{^(F r )} (44) 

exists, is positive and finite, where 

i>(Y) = (2 " f) ^(0) + + f e*<°> + w(f (0)) + ^ Y-(-t)) dt. (45) 

We will drive the diffusions (|42|) with the same Brownian motion B(t). Then for T\ > T 2 
we have Y Ti (t) > Y t . 2 (t) for r £ [T 2 , 00) as this holds for r = — T 2 and the domination is 
preserved by the evolution. 

We also consider a nonnegative-valued diffusion Z{t) given by the SDE 

dZ = r{Z)dt + dB 

which is reflected at and whose drift term is equal to 

r(v) = -Tz? + c x . (46) 
lb 



By Lemma[7]the term h (y, T + r) in (|43)> is bounded if — T < r < 0. Thus we can choose 
the constant c\ in (|46|) so that 

r(z) > sup h(r,y). (47) 

{r<0,0<?/<z} 
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Since Z and Y are driven by the same Brownian motion, if Z, Y > then Z — Y evolves 
according to 

d(Z-Y) = [r(Z)-f(t,Y)]dt. 

By (|47|) this means that if Z(r ) > Y(t ) for a r < then this ordering is preserved by the 
coupling until time 0. 

Consider the process Z in its stationary distribution. Then Z(—T) > Yt(—T) there- 
fore Z dominates Yr on [— T, 0]. For every fixed r < the random variables Y t (t) are 
increasing in T and bounded by Z{r) so 

Y^t) = lim Y t (t) 

T — >oo 

exists and is dominated by Z(r). The function pi(x) is continuous sopi(Y T ) — > pi(Yoo). By 
((45]) we have 

$(Y T ) = a(Y T (0)) + / <p(t, Y(-t)) dt 



where a is continuous and y) can be dominated by a function <j)(y) which has a finite 
integral in [0, oo). Hence ip(Y T ) — > ^(Y^) and 

q T = e^ Pl (Y T ) ^ goo = e^Vi(^co), as T ^ oo. 

Using Lemma [5] to estimate we get 

q T < ce xp|(2-/3/2)y T (0) + + ^°)-^e^ )} < c' X (^t(0)), 

where x(y) — ex P {(f — §[) eV }- If we prove that Ex(Z(0)) < oo then the dominated 
convergence theorem will imply 

Eq T -> Eg^ < oo, (48) 

and the existence of the limiting constant Kp will be established. 
The generator of the reflected diffusion Z is given by 

Cf= 1 -f" + fr, 



for functions / defined on [0, oo) with f(0+) = 0. (|Revuz and Yori (|1999l) . Chapter VII. §3.) 
Partial integration shows that if (logg)' = 2r and /'(0+) = then f °° Cf(x) g(x)dx = 
which means that 

g{z) = cexp(-(3/8e z + 2c x z) 
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gives a stationary density. Since J °° x( z )g( z ) dz = E%(Z(0)) < oo, the convergence ((48]) 
follows. This shows that 

Kfi = Egoc = E [e^pi^O))] < oo. 
The only thing left to prove is that k@ = E is not zero. The definitions of q and t/> yield 

goo>c Pl (y oo (0)) e ^/ 2 ) y -( ) + . 

By Lemma|5]the function pi(-) is positive. Since i^o(0) is a.s. finite and we get that Eg^ > 
which completes the proof of Theorem [U □ 
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